Pressure Fronts in ID Damped Nonlinear Lattices 
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The propagation of pressure fronts (impact solutions) in ID chains of atoms coupled by anhar- 
monic potentials between nearest neighbor and submitted to damping forces preserving uniform 
motion, is investigated. Travelling fronts between two regions at different uniform pressures are 
found numerically and well approximate analytically. It is proven that there are three analytical 
relations between the impact velocity, the compression, the front velocity and the energy dissipation 
which only depend on the coupling potential and are independent of the damping. Such travelling 
front solutions cannot exist without damping. 

I. THE MODEL 

This study was motivated for understanding sonoluminescence (see refs. ini) observed for example when water 
saturated with rare gas is submitted to an intense ultrasonic field. Spherical bubbles of rare gas expand and collapse 
periodically at supersonic or nearly supersonic speed emitting simultaneously a short and intense broadband light pulse 
at impact on the bubble core (sonoluminescence). We suggested more generally that light emission is systematically 
generated at strong enough impacts still below the range for generating plasmas^ but providing they become highly 
nonlinear. This situation occurs when the hard core of the atoms (which prevent volume elements to become negative), 
is involved that is for supersonic or nearly supersonic impacts. We suggest now a new remark for amending our early 
theory. Since condensed matter is made of bonded charged particles, sharp accelerations of these charges at strong 
enough impact may be sufficient to generate an intense (Abrahams-Lorentz) electromagnetic (em) radiation visible 
as sonoluminescence and generating damping (note that a charged hard sphere model would produce a diverging 
radiation!). We do not discuss here the physical validity of this suggestion but we only focus on some preliminary 
mathematical aspects of shocks in simple ID model with damping. 

FPU lattices are ID chains of atoms coupled by anharmonic springs with Hamiltonian 

where u„(i) is the scalar coordinate of atom n, pn ~ u„ is the associated conjugate variable and V{v) is the coupling 
potential which depends on the distance w„ = ii„_(-i — u„ between nearest neighbor atoms n + 1 and n. 

We are interested in moving pressure fronts between two regions at different pressure obtained from initial conditions 
corresponding to an impact which are for example: Un(0) = nv and u„(0) = — Vp for n > 0, Un{0) = +Vp for n < 
and vo(0) = 0. The chain at equilibrium is initially moving uniformly for positive n and at the opposite velocity for 
negative n (Vp is the impact velocity). Since by symmetry arguments, atom remains immobile at all time (uo(t) = 0), 
it is equivalent to consider that the positive part of the chain n > 0) impacts a fixed rigid wall at uo(0) = or in 
the framework of the center of mass of the half chain that atom in the chain initially at rest is pushed by a piston 
uo{t) = Vpt. Studies were already devoted to pressure fronts in lattice models without damping^ii^. The same model 
with damping preserving uniform motion (as in^ for a Toda potential) is described by eqs. 

Un - 7(u„+i + ?i„_i - 2u„) - V'(u„+i - u„) + V'(u„ - u„-i) = or (2) 

Vn -jA.Vn - A.V'iVn) = (3) 

where operator A is defined by A.F„ = Fn+i + Fn-i — Fn and 7 > is the damping constant. A solution of eq[3] 
is Vn{t) = V where v is an arbitrary constant . The linearized equations for |w„(t) — v\ small at 7 0, yield plane 
wave solutions Vn{t) = v + Acos{qn — uj{q)t — a) where uj'^iq) = 4V"(v) sin^ q/2 « s'^{v)q^ for small q. For avoiding 
more complex situations where the front breaks into several fronts ( phase separation), it is convenient (and physically 
reasonable) to assume that the square of the sound velocity s^iv) = V"(u) is a monotone decreasing function of v 
that is V"{v) < 0. 



2 



II. CONTINUOUS MODEL 



Since continuous model are often used for describing fluids, we investigate first the continuous version of this model 
which will appear physically inconsistent but nevertheless will reveal interesting features. Assuming — u„| small 

that is u„ {t) is a slowly varying function of n = x (while the variation of is not necessarily small), the PDE 



d^u dV'{—) 

~^dxm - ^ ° °q^-^l^-tly (4) 



describes u„(i) = u(x,t) or Vn{t) = v{x,t) = Eq[5] exhibits exact step front solutions v{x,t) ~ f-oo + (v+oo — 
?;_oo)Y(x — ct) (with fast variation!) corresponding to u{x, t) = v^^oX + (w+oo — V-oo){x — ct)Y(x — ct). (Y(.t) is the 
Heavyside function (Y(x) = for a; < and Y(x) = 1 for a; > 0). The square of the front velocity c 



2, ^ V^(^'+oo)- V'(^-oo) ... 
C (u_cx),W+oo) = (6) 

W+oc - V-oo 

is a only function of the atomic compression at infinities w+oo and v-oc- The difference Vp ~ ?i(— oo, t) — u(+oo, t) is 
the impact velocity 

Vp = c(u+oo - W-oo) (7) 

Considering the energy of a finite but long part of the chain containing the front $(f) = {^ii? + V{^)) dx , 
we readily obtain for any solution u{x,t) of eql4]with smooth second derivatives that the rate of energy variation of 
the system per unit time $ = $o where $0 = V {u{L , t))u{L , t) — V'{u{—L,t)u{—L,t) is the power delivered by the 
pressure at the edge of the system (energy conservation). When the second derivatives of u(x, t) are not smooth but 
involve Dirac functions, there is generally no energy conservation! An explicit calculation of the energy $ for the step 
front solution (using cq. ^) yields the dissipated power 

D = ^-^0 = -C (v{v+oo) - V(W-oo) - ^{v+oo - l^-oo)(V'(«+oo) + V'(z;_oo))) (8) 

which may be written again in a simpler form as 

D = - V-oo)^S{v+oo,V-oo) (9) 

where S{x,y) = ^^q^:^^ x-derivative of R{x,y) defined by the equation V{x) = V{y) + (x ~ y)V'{y) + ^{x — 

?/)2 (V"(y) + Rix, y)). We remark that 




is non vanishing only when potential V is anharmonic. Then, the above assumption V"'{x) negative readily implies 
{x — y)R{x, y) < and (a; — y)S{x, y) < 0. Thus, we have spontaneous energy dissipation when c{v+oc — Voo) > or 
energy creation in the opposite case! 

Because of this lack of energy conservation, this continuous model is not physically acceptable for a correct de- 
scription of impacts. Where this energy would go (or come from)? This continuous equation at zero damping has 
another flaw we do not detail here. There are smooth initial conditions which spontaneously develop a singularity (a 
divergence of the derivative ^ (x, t)) within a finite time by the standard Rankine-Hugoniot mechanism. Actually, the 
pressure fronts in lattice models without damping are not stcplike but exhibit a puzzling behavior with an expanding 
intermediate region and backward oscillations extending to —00 where nearest neighbor atoms are in antiphase^. 
Despite the inconsistencies of this continuous model, we prove now that formula [71 [6] and [8] which relates the impact 
velocity Vp = — u+00 + u-oo, the compression w+00 and V-00 at infinity, the front velocity c and the power dissipated 
by the front, are the same for the discrete case when damping is present but independently of this damping. 
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III. PRESSURE FRONTS IN DISCRETE MODEL: EXACT RELATIONS 

We assume in the discrete model [2] with 7 > 0, there are traveUing front solution Vn{t) ^ g{n — ct) described by a 
smooth stcplike hull function g{x) fulfilling \mi.^^±oo g{^) — '^±00 , ^^''^x~^±oo g'i^) — and 

c^g" + -fcA.g'ix) - A.V'(g(x)) (10) 

(Operator A is defined by A.F{x) = F{x + 1) + F{x - 1) - 2F{x)). 

For proving eqlll we consider the length L-N,+Nit) = UN{t) — u^N{t) — '^^jsr^ Vn{t) = X]-Ar^.9("- ^ ^t) of 
the chain between far sites —N and +N. We have liraN^+oo L-N,+N(t) = liniAr^+oo (u+w — u-n) = u+00 — 
■it„oo = —Vp. Otherwise, because of the existence of the hull function g{x), the rate of variation of the length 

c(i_Ar,+Ar(i + 1/c) — L_jv^+7v(i)) between time t and t + 1/c is also equal to c{—g{N — ct) + g{—N — 1 — ct)) = 

c(— ujv + v-N-i)- Then, limjv^+oo ^i-Af,+Af(i)^ = —c{v+oo — ^'-oo) = —Vp which yields eqlT] 

For proving eq. [SI we consider the momentum A/_7v.+Af(i) = X)-jv "^nit) of the chain between far sites —N and +N. 

EqUyields M_7V,+Af = j{vN-v-N-i) + y{vN)-V'{v-N-i) and then M = liniN^+oo M_n,+n = V'(u+oo)- V'(v_oo)- 
Eql2]also yields it„ — G{n — ct) where G{x) = —cj{g'{x) — g'{x — 1)) + V'{g{x)) — V'{g{x — 1)). Time integration 
yields = F{n — ct) + a„ where F{x) ~ -~G'{x)/c and a„ is time-constant. Since Vn = —cg'{n ~ ct) ~ iin+i — Un = 
F{n + 1 — c^) — F(n — ct) + a„+i — a„, it comes that a„+i — a„ which is time constant, should be only a function 
of 71 — ct. It is thus also independent of n. Then it comes out a„+i — a„ = when considering the limits n —>■ ±00. 
Thus, the rate of variation over the period of time 1/c of the momentum c(Af_jv,+Ar(t + 1/c) — M_jv,+Af(t) ~ 
cE^jv(^(" - 1 - ct) - F{n - ct)) = c{-F{N - ct) + F{-N - 1 - ct)) = c{-u+n + ii-N-i) becomes equal to cVp 
when N +00. It should be equal to M = limjy^+oo l^-t-N,+N = ^'('f+oo) ^ V(w_oo) which is time constant when 
N +00. Combining this equality with ([7]) yields eqlHl 

Finally, for proving eq[51 we consider the energy of the chain $_jv,+Ar — '^n=-N h^n~^^^N^ V(u„+i — u„) between 
sites -N and +N. We readily obtain <i>_Ar,+Ar = Yln=-N ~ V'(w„+i - u„) + V'{un - Un-i)) Un + V'{un+i - 
un)un — V(u_Ar — U-N-i)u~N ■ We define <i>^]^ _).jy = V'{un+i — un)un — V'{u-n — u^n^i)u^n as the power 
provided by the external pressure to the finite chain. We have for the infinite chain <i>'-'^' = limAr^-|_oo +jv ^ 

V'(i;+oo)u+oo - V'('y-oo)u-oo- Then, using cqEl we readily obtain ^^n,+n - *^lr,+Ar = ~7 1]n=-jv("n+i ~ "n)^ + 
(mtv+i — un)un — (u-N — u_Ar_i)u_jv Sincc limyv^+oo u±N = 0, we obtain that the power dissipated by the damping 
force in the infinite system D = liniN^^ooi^ -n,+n — +n) < is 

+00 

i?--7 5I(""+i-"")' (11) 

—00 

This dissipate power may be calculated differently sincc wc assume that the solution is described by a hull function 

Vn = g{n - ct) implying u„ = F{n - ct). Then, we have $_jv,+Af(t) = Y.n=-N ^F^in - ct) + Y.n=-N ^i9{n - ct)) 
so that we readily obtained similarly as above, the rate of variation over an interval of time 1/c of the energy 
c{^-N,+Nit + 1/c) - ^-N.+Nit)) = c{\{F^{-N - 1 - ct) - F^{N - ct)) + V{g{-N - 1 - ct)) - V{g{N - 1 - ct)) which 
for N —f +00 becomes c{^{u^_^ — u\^) + V(w-oo) — V[v+oo)) = $ = hmAr_^+oo *i'-Af,+JV which is time constant. 
Consequently, we obtain for the infinite chain, D = $ — d?'*^^ = c(^(u_oo ~ u+oo){u-oa + ^+00) + V(w_oo) — V(f+oo)) — 
V'(v+oo)w+oo + V(w_oo)w_oo- Using Vp = U-oo ^ and eqsl7|and[ni the dissipated power D defined by eafTTl is 
found to fulfill eqlSl 

As a consequence, there is no stationary travelling front solutions g{x) in the harmonic system with damping 7 7^ 
because eq|8] yields D = while eq[Tl] yields D ^ 0. There is also no stationary travelling front solutions when the 
system is anharmonic and 7 = because eq|8] yields D ^ while eq[Tl] yields D = 0. Otherwise, since D has to be 
negative, front solutions may only exist when they propagate from larger toward the lower pressure region (c > 
when u_oo < ^'+00)- 

Assuming the existence of a hull function, g{x), it is straightforward to extend the same proofs for formula [71 [5] and 
[5] to any other kind of damping forces preserving the uniform motion of the chain for example the Abrahams-Lorentz 
force proportional to the third time derivative m„. 




IV. NUMERICAL CALCULATION OF PRESSURE FRONTS 



The hull function g{x) at nonvanishing damping can be numerically calculated at computer accuracy. We start 
form initial conditions where a half chain is at rest u„(0) = nv^ao and u„(0) = for n > while the edge site 
uoit) = Vpt is constrained to have a uniform motion at velocity Vp > (piston or impact velocity). The front velocity 
c{Vp, u+oo is determined as a function of the impact velocity Vp and u+oo through cqs [6] and [7] which yields the implicit 
equation c = (V"(u+oo) — ^'(w+oo — Vp/cj) /Vp. It is a monotone increasing function of Vp because V" is assumed to 
be monotone decreasing. At Vp = 0, c^(0, w+oo) = V"(u+oo) = s^(w+oo becomes the sound square velocity. The map 

/ K+i(i/c)} \ / K(o)} \ 

V {Vn+l{l/c)} [ {Vnm ) 

defined by integration of eql2] or [3] over the period of time 1/c and a shift of the indices, is iterated by numerically 
from the impact initial conditions X defined above. It is found that when the damping constant 7 > is non zero, 
TP{X) systematically converges for p +00 to a fixed point which corresponds to a solution Vn(t) = g{n — ct). Plots 
of an example of calculation of this hull function is shown fig[l] for several damping constants. At strong damping 
7 ~> 1, convergence is obtained within few iterations only and the hull function is step like. For smaller damping, 
the rate of convergence slows down while the stationary front solution develops backward oscillations which diverges 
at 7 = 0. At zero damping, there is no convergence at all. Actually, this problem was already investigated in the 
literature^. 

At nonvanishing damping, g{x) may be well approximate as a solution of a differential equation describing an 
anharmonic oscillator with damping. We apply operator 1 — Q where Q = Qi j^) with Q{z) = 1 — 2(cosh2-i) ~ 

^' il>o^'"/S"ti))i ^ Ep>o'?"^'" = 15(1 - 15^^ + ■•••) to the left member of cqM Since A ^ L{^) where 
L{z) = 2(coshz — 1) = jTi^py, eqfTOl becomes ^ {c^{l — Q).g + ^cg' -~V'{g)) = and after two integrations 
^ 7C.9' + — c?g ~ ax + b where a and b are two arbitrary constants. Since we search for (physical) solutions 
g{x) which arc bounded at ±00, wc must have a = 0. Then, defining potential W(g) = V{g) — \c^g^ — bg, eg [TUl takes 
the form 

c2Q..g-7cg' + W'(.g) = or (12) 
^5" - 7C5' + W'(<?) « (13) 

for the lowest order approximation Q{z) « j^^^- Eq llSl mav be viewed as the equation of a negatively damped particle 
with coefficient —7c and mass c^/12 in the effective potential W(g). There are non diverging solutions for x ±00 
only when }V{g) has at least two extrema and then the solution is asymptote to each of these extrema. Actually since 
W'{g) = V"{g) — (? is monotone decreasing, W((7) has at most a maximum at g = w+oo > ''^-oo and a minimum at 
g — w_oc which are determined by b and (? and then eq[5]is fulfilled. FigU] (right) shows a fit of the hull function 
obtained as the separatrix solution of this equation 1131 such that limj;^_oo .9(2^) = '^-oo is the minimum of Wix) and 
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lim2;^_i_oo gix) = f+cx; is the maximum. The error is negUgible at large 7 but appears mostly as a phase shift in the 
tail for small 7 visible fig[T] right. 

In summary, the most important result of this paper is that the energy dissipated at an impact is independent of 
the physical origin of the damping and of its value providing it preserves the translational motion of the system (as 
it should in physical models). It only depends through formula ([9]) on the anharmonic part of the potential between 
the compressions ahead and backward the front and is proportional to the front velocity. This result would formally 
determine the emitted power of sonolumincscence if one believes our physical interpretation. It is then straightforward 
to check that the emitted power is negligible for nearly harmonic impacts but in principle could approach 100% of 
the input power when V is a hardcore potential diverging at some Vc < v and for strong impacts where u_oo ~ ^c- 
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